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Abstract: We report the results of a successful 12 hour 22 GHz VLBI experiment using a heteroge- 
neous network that includes radio telescopes of the Long Baseline Array (LBA) in Australia and several 
VLBI stations that regularly observe in geodetic VLBI campaigns. We have determined positions of 
three VLBI stations, ATCA-104, CEDUNA and MOPRA, with an accuracy of 4-30 mm using a novel 
technique of data analysis. These stations have never before participated in geodetic experiments. We 
observed 105 radio sources, and amongst them 5 objects which have not previously been observed with 
VLBI. We have determined positions of these new sources with the accuracy of 2-5 mas. We make 
conclusion that the LBA network is capable of conducting absolute astrometry VLBI surveys with 
accuracy better than 5 mas. 
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1 Introduction 

The method of very long baseline interferometry (VLBI) 
first proposed by Matveenko, Kardashev, & Sholomit- 
sky (1965) is widely used for geodesy and absolute as- 
trometry. The first dedicated geodetic experiment on 
1969 January 11 provided results with 1 meter accu- 
racy (Hinteregger et al. 1972). In the following decades 
VLBI technology has flourished, sensitivities and ac- 
curacies have improved by several orders of magni- 
tude and arrays of dedicated antennas have been built. 
Nowadays, a network of ~30 stations spread over the 
globe participates regularly in observing programs for 
the Earth orientation parameters (EOP) monitoring, 
determination of source coordinates, deriving site po- 
sitions, and monitoring their changes. These activities 
are coordinated by the International VLBI Services for 
Astrometry and Geodesy (IVS) (Schliiter & Behrend 
2007). 

However, the station distribution over the globe is 
non-uniform. The majority of radio telescopes are lo- 
cated in the northern hemisphere. This non- uniformity 
results in a disparity in the density of source cata- 
logues: the number of sources with positions deter- 
mined at a milliarcsecond level of accuracy in the dec- 
lination zone [—90°, —50°] is a factor of 5 less than in 
the zone [+50°, +90°]. Disparity in station distribu- 
tion also results in the appearance of systematic errors 
in the EOP time series derived from analysis of VLBI 
observations at this network. Therefore, an increase 
of the number of observing stations in the southern 
hemisphere that are able to participate in VLBI ob- 



serving campaigns under geodesy, absolute astrometry, 
and space navigation programs is very important. 

There are seven radio telescopes in Australia that 
have VLBI recording equipment: HOBART26, parkes, 
tidbin64, DSS45, Australia Telescope Compact Array 
(atca), CEDUNA, MOPRA, and they potentially can 
participate in VLBI observing campaigns for astrome- 
try and geodesy. The first four stations are equipped 
with Mark-4/Mark-5 data acquisition system and they 
have participated in many IVS observing programs. 
The last three stations have the Long Baseline Array 
Data Recorder (LBADR) data acquisition system that 
is not directly compatible with Mark-4/Mark-5 and 
they have never before participated in experiments for 
geodesy and absolute astrometry. Using all stations 
together is a challenging problem. 

The possibility to use sensitive antennas CEDUNA, 
ATCA, MOPRA for space navigation and astrometry 
would significantly boost our capabilities to observe 
objects in the southern hemisphere. On 2005 Jan- 
uary 16 stations ATCA-104, CEDUNA, MOPRA and 
parkes joined the global VLBI network that also in- 
cluded GBT, eight VLBA stations kashima, seshan25, 
and URUMQI in the Huygens probe tracking observa- 
tions during its descent in the atmosphere of Titan 
(Pogrebenko et al. 2004; Witasse et al. 2006). De- 
tection of the Huygens probe signal on the North 
America-Australia baselines was crucial for the recon- 
struction of the probe trajectory with ~ 1 km accuracy 
at Saturn's distance of 1.3 billion km, although the 
uncertainty of the LBA station a priori coordinates in- 
troduced systematic errors in the probe's positioning. 
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A better a posteriori determination of the site coor- 
dinates was required to improve the accuracy of the 
probe's trajectory reconstruction. 

The Long Baseline Array incorporates 5 main an- 
tenna: the parkes 64 m, the MOPRA 22 m, the HO- 
BART 26 m, the CEDUNA 30 m, and the Australia Tele- 
scope Compact Array which consists of six 22 m anten- 
nas. When ATCA participate in VLBI observations, 
it can run as a phased array or, alternatively, only 
one of its antennas can be used. Antennas from the 
NASA's Deep Space Network at the Canberra Deep 
Space Communication Complex, tidbin64 and DSS45, 
regularly join the array as does the hartrao 26 m dish 
in South Africa. The array is used for source imaging, 
differential astrometry and other applications. 

Imaging weak target sources and precise differen- 
tial astrometry require a catalogue of compact calibra- 
tors with the density at least one object within a circle 
of 2-3° of any direction. A catalogue of such calibra- 
tors was created from analysis of the VLBA Calibra- 
tor Survey observing campaign (Beasley et al. 2002; 
Fomalont et al. 2003; Petrov et al. 2005, 2006, 2008a; 
Kovalev et al. 2007). However, the VLBA cannot see 
objects with declinations below —50°. The extension 
of the Calibrator Survey to the southern hemisphere 
should be derived from analysis of dedicated LBA ob- 
serving campaigns. 

For investigating the possibility of using the LBA 
antennas for precise geodesy and absolute astrometry 
we ran two observing sessions: a 1.5 hour long test 
VLBI experiment and a 12 hour geodetic VLBI exper- 
iment. The goal of these experiments was 1) to de- 
termine positions of ATCA, CEDUNA, MOPRA with the 
position accuracy 1-5 cm; 2) to evaluate the feasibility 
of running absolute astrometry observing campaigns 
with the heterogeneous LBA network. 

In this paper we present results of these exper- 
iments. In Section 2 we describe the experiments, 
scheduling strategies and the hardware configuration. 
In Section 3 we describe the way how the hetero- 
geneous data were transformed to the same format, 
transmitted to the correlator and correlated. In Sec- 
tion 4 we describe the algorithm for data analysis and 
present our results. Concluding remarks are given in 
section 5. 



2 Observations 

On 2006 May 16, stations atca-104, ceduna, mopra, 
HOBART26, parkes ran a 1.5 hour long fringe test ex- 
periment at 22 GHz. The purpose of the fringe exper- 
iment was to test the data path, to test the frequency 
selection, to detect fringes, and to test data analy- 
sis procedure. Stations atca-104, ceduna, mopra 
recorded with the LBADR system. Data from these 
stations were converted into Mark-4 format using soft- 
ware developed at Metsahovi Radio Observatory for 
the Huygens VLBI experiment, then the data were 
test-correlated at the Joint Institute of VLBI in Eu- 
rope, Dwingeloo (JIVE) in Dwingeloo, the Nether- 
lands, and finally correlated at the Mark-4 correlator 
at the MPIfR in Bonn, Germany. Fringes were found 



at all baselines, except baselines with ceduna. Re- 
sults from this test run allowed us to determine coarse 
positions of ATCA-104 and MOPRA, which helped to 
improve noticeably the Huygens VLBI trajectory re- 
construction and encouraged us to run to a full-scale 
22 GHz LBA experiment in the geodetic mode. 

The telescopes atca-104, ceduna, hartrao, ho- 
BART26, MOPRA, parkes, and SESHAN25 (Figure 1) 
took part in the 12 hour experiment on 2007 June 24. 
Only one antenna of the 6-element ATCA array par- 
ticipated in the experiment. Antennas of the ATCA 
array can move along railroad tracks and can be po- 
sitioned at one of 45 fixed station pads within 6 km 
range. The ATCA antenna that was used in the obser- 
vations was positioned at the station pad CATW104, 
and it is referred to as ATCA-104. 

This experiment was very different from routine 
geodetic VLBI experiments. 

• We chose to observe at the K band (22 GHz). 
Traditionally, geodetic experiments are made in 
a dual-band mode when emission at S band 
(2.3 GHz) and X band (8.4 GHz) is recorded 
simultaneously. The linear combination of de- 
lays at the X and S bands is almost entirely free 
from the contribution of the ionosphere. The an- 
tennas ATCA-104, CEDUNA, MOPRA do not have 
dual-band S/X receivers, so we are limited to 
only one band. The highest frequency band 
that all antennas can receive is 22 GHz. Since 
the ionosphere contribution to group delay is 
reciprocal to the square of frequency, the iono- 
sphere contribution at the K band is one order of 
magnitude smaller than the contribution at the 
X band. Our analysis of K band astrometry ob- 
servations with the VLBA and several K band 
geodetic experiments with the Japanese VLBI 
interferometer network VERA showed that the 
ionospheric contribution at this frequency dur- 
ing the period of low solar activity is negligible, 
and the overall quality of geodetic results turned 
out even better than of dual band S/X observa- 
tions. 

• The digital backend used at ATCA, MOPRA and 
PARKES allows recording of only two pairs of 
adjacent 16 MHz intermediate frequency chan- 
nels (IF) in right circular polarization and the 
two pairs of adjacent frequency channels in the 
left circular polarization, ceduna can record 
only two pairs of adjacent IFs in one polar- 
ization. Although the IFs are recorded inde- 
pendently, the system effectively records emis- 
sion from four channels and always records dual 
polarization (except ceduna), effectively giving 
only two 32 MHz channels at different frequen- 
cies at each polarization. This differs signifi- 
cantly from other VLBI backends such as the 
Mark-4/Mark-5 and K5 systems which support 
up to 32 frequency channels and the VLBA 
recording system which supports up to 8 fre- 
quency channels. Traditionally, from 4 to 16 
channels per band are allocated. The frequency 
sequence is selected to minimize the power of 
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Figure 1 The network of VLBI stations that participated in the experiment on 2007 June 24. Stations with coordinates 
known from analysis of prior experiments are shown with big disks. 



sidelobes at the delay resolution function and 
for avoiding possible radio interferences. No ge- 
odetic experiments have ever been made with 
only two channels. 

The frequency sequence that we have chosen for 
this experiment is shown in Table 1. The fre- 
quency difference between two effective 32 MHz 
channels is 224 MHz. Since the accuracy of 
group delay determination observed at two ef- 
fective intermediate frequency channels is recip- 
rocal to the frequency separation, the wider the 
frequency separation, the better. From the other 
hand, the wider the separation of two effective 
IFs, the smaller group delay ambiguity spacing 
that is also reciprocal to the minimum frequency 
separation between separate frequency channels. 
The group delay ambiguity spacing in tradi- 
tional VLBI experiments ranges from 28.5 ns 
to 200 ns. The spacings were 1/224 MHz = 
4.5 ns in our experiment! Although the theoret- 
ical model for VLBI path delay has truncation 
errors less than 0.1 ps, the accuracy of delay 
prediction is limited to the accuracy of a priori 
parameters of the model. The major uncertainty 
of the a priori model arises from the contribu- 



tion of the wet component of the tropospheric 
path delay that currently has a typical accuracy 
of 3-10 ns — higher than the group group de- 
lay ambiguity spacings. Therefore, traditional 
methods of data analysis that use group delays 
would fail. 

Table 1 The range of sky frequencies at four frequency 
channels in the observing session of 2007.06.24, in GHz. 



FC1 


22.300 - 


22.316 


FC2 


22.316 - 


22.332 


FC3 


22.524 - 


22.540 


FC4 


22.540 - 


22.556 



• All stations, except CEDUNA, recorded both right 
circular polarization (RCP) and left circular po- 
larization (LCP). Traditionally, only the RCP 
is recorded in geodetic experiments. Since the 
thermal noise of RCP and LCP data is indepen- 
dent, the use of LCP data reduces random errors 
of parameter estimates. 

• The experiment had a heterogeneous setup. 
Data were recorded at 512 Mbit/s, 2 bit sam- 
pling, 16 MHz bandwidth per IF, 8 IFs, dual 
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circular polarization (except from CEDUNA that 
can record only RCP). HOBART26, hartrao and 
Seshan25 used a Mark-4 data acquisition termi- 
nals with the setup shown above and wrote data 
onto disks using the Mark- 5a VLBI recorder 
(Whitney 2005). The rest of the antennas used 
digital data acquisition system and the LBADR 
VLBI recorder 1 which is a derivative of the 
PCEVN 2 (Parsley et al. 2003). Data were re- 
corded as 16 parallel bit streams. The LBADR 
uses different physical media for recording than 
the Mark-5a. 

The LBADR system uses the VLBI Standard In- 
terface (VSI-H) (A. Whitney 2000, unpublished inter- 
nal document 3 ) signaling and connector conventions. 
Output bit streams from the VLBI digital backend are 
converted to the VSI-H using the universal VSI con- 
verter board called "VSIC" and captured into regu- 
lar files on a server-grade PC run under Linux with 
a VSI-H-to-PCI record board called "VSIB". To en- 
sure sufficiently high data record rates, Linux software 
"raidO" disk sets were used to store the recorded files. 
This is different from Mark-5 which uses proprietary 
disk modules as its physical media for recording. 

To achieve compatibility with the correlator Mark- 
5 playback systems, two software solutions were used. 
The LBADR data of the first 1.5 hour fringe test run 
were software transformed to Mark-4 format, main- 
taining timecode information, using custom software 
running on a computer at JIVE. The resulting files 
were transferred to a Mark-5a unit to record disk mod- 
ules for the Bonn correlator. 

For the subsequent 12 hour run, software was de- 
veloped that converts on-the-fiy the raw LBADR data 
to Mark-5b format then sends the data using a tuned 
TCP stack over the network directly to files residing on 
a data server at JIVE. Mark-5b compatible disk mod- 
ules were created and shipped from JIVE to Bonn for 
correlation. For future such experiments, the LBADR 
system will write Mark-5b data format directly, but 
will still not write to Mark-5 diskpacks. 

2.1 Scheduling 

The input list for automatic scheduling consists of 
129 compact radio sources. This list was generated 
by merging two lists for different declination bands. 
Sources in the declination band [—30°, +40°] were se- 
lected from the catalogue of the 252 objects observed 
at K and Q band with the VLBA in 2002-2006 in 
the framework of the K/Q astrometric survey (Fey 
et al. 2005; Jacobs et al. 2005). They have the me- 
dian K band correlated flux density > 0.4 Jy at base- 
lines longer than 5000 km. Sources in the declination 
band [—90°, —30°] were taken from the list of sources 
that were a) observed in geodetic VLBI experiments, 
b) observed and successfully imaged in astronomy ex- 
periments (Ojha et al. 2004), and c) have the median 
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X band correlated flux density > 0.4 Jy at baselines 
longer than 5000 km. 

Automatic scheduling was made with the program 
SKED developed at the NASA GSFC by N. Vanden- 
berg and J. Gipson (Petrov et al. 2008b). The algo- 
rithm for automatic scheduling generated a sequence 
of scans, i.e. intervals of time of 120 s duration when 
all antennas of the array, or a portion of the array, are 
pointed to a specific source and record the emission. 
Since the list of sources has 129 objects, there exists 
a very large number of possible sequences. SKED se- 
lected the sequence that optimizes the sky coverage at 
each individual station and minimizes the slew time 
needed for re-pointing antennas at the next source of 
the sequence. At each step of generating the sequence 
of scans, the software finds the distance of a candidate 
source from all previous sources scheduled within one 
hour. The candidate source with the largest minimum 
distance from all previous sources gets the highest sky- 
coverage score. SKED computes the slewing time and 
assigns weights to all candidate sources according to 
their scores based on sky-coverage, slewing time and 
some other optimization criteria. The source with the 
highest weight is put into the schedule. The scheduling 
algorithm adjusts weights to low elevation sources in 
such a way that each station observes at least one low 
elevation source every hour. The elevation cutoff was 
12° for all the stations, except parkes, that cannot 
slew lower than 31° above the horizon. 

SKED automatically selected 100 sources. In ad- 
dition to them, five bright flat-spectrum sources from 
the Parkes Quarter- Jansky catalogue (Jackson et al. 
2002) were inserted in the schedule manually at sta- 
tions ATCA-104, CEDUNA, MOPRA, HOBART25, PARKES, 

two scans each. These sources have declination below 
— 55° and have never been before observed with VLBI. 
The goal of including these sources in the schedule is 
to evaluate the feasibility of using the LBA for a search 
of new compact radio sources that can be used as cal- 
ibrators for phase referencing observations. 

3 Data Transmission and Cor- 
relation 

Because of the media incompatibility of the LBADR 
and the Mark-5 format, as well as for convenience, data 
were transmitted using high speed networks between 
Australia and Europe. The data transmission was per- 
formed directly from the recorder PCs at the observa- 
tories in Australia, to a PC located at the JIVE. This 
was to utilize a dedicated 1 Gbps network which had 
previously been set up for eVLBI demonstrations as 
a part of the Express Production Real-time e-VLBI 
Service (EXPReS). The transmission was made using 
custom software which uses the TCP network proto- 
col and has been tuned to efficiently use long, wide 
bandwidth networks. This software also performed an 
on-the-fly conversion of the LBADR data format to 
Mark-5b. This simply requires removing the LBADR 
headers and replacing them with Mark-5b headers (re- 
taining time stamps) as the baseband bit format can 
be processed by Mark-4 correlators. 
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After transmission, the Mark-5b data were copied 
from the PC at JIVE to a Mark-5 disk-pack and ship- 
ped to Bonn correlator. Fringe checking was run on 
the EVN Mark-4 data processor at JIVE before the 
disk-packs were shipped to Bonn for final correlation. 
Those stations, that recorded with the Mark-5a data 
acquisition system, shipped the diskpacks with data to 
the correlator directly using air mail. 



3.1 Correlation 

The experiment was correlated at Bonn since the Bonn 
correlator has previously been thoroughly tested for 
astrometric and geodetic applications. 

The Mark-4 correlator (Whitney et al. 2004) was 
configured to produce 512 lags and a 0.5 s accumu- 
lation period, to have as wide as possible delay and 
delay-rate fringe search window to allow for possible 
large station coordinate errors and source position er- 
rors. For comparison, the routine geodetic experi- 
ments use only 32 lags and accumulation periods of 
between 2 and 4 s since the errors in the a priori station 
and source positions are small. In this configuration, 
the correlator computing capacity is sufficient to pro- 
cess only four stations. Thus, we correlated six passes 
to form all the baselines between the seven antennas 
(with some redundancy). The selection of stations for 
correlation in each pass was guided only by the num- 
ber of available Mark-5a and Mark-5b units at Bonn, 
but paying attention not to transpose data from the 
same antenna when correlating redundant baselines. 

Since the LBA stations were not controlled by the 
Mark-4 field system software and therefore, they did 
not produce the observation log file required for corre- 
lation. Instead, the information required for the cor- 
relator control files, such as the scan name and scan 
times were recovered by comparing information from 
the Mark-5 module to the schedule file. 

A trial correlation was performed to check the cor- 
rectness of the correlator control files, for example 
clock offsets, frequencies, polarizations and track as- 
signments. Quality control checks after each pass were 
made using the Haystack observatory post-processing 
system (HOPS) software package. 

The stations did not inject phase calibration tones, 
so phase self calibration was applied to all stations in 
order to remove IF-dependent phase offsets at those 
stations with analogue data acquisition terminal (ho- 
BART26, HARTRAO, and SESHAN25). In the 17% of 
cases where the source was not detected on a base- 
line, those scans were re-correlated to ensure that these 
non-detections were not correlator artifacts. 

Finally, the multi-band group delays and delay 
rates were adjusted using the fringe fitting algorithm 
implemented in HOPS program Fourfit in such a way, 
that applying these corrections to group delay and de- 
lay rate would maximized fringe amplitudes over each 
scan at each baseline independently. 



4 Data Analysis 

The fringe fit algorithm provides the estimates of 
fringe amplitudes, fringe phases, phase delay rates, 
narrow-band group delays, and multi-band group de- 
lays. Analysis of amplitude information provides the 
estimates of the correlated flux densities of observed 
sources and the sensitivity of antennas. Analysis of 
fringe phases and delays allows us to estimate station 
positions and source coordinates. 

4.1 Position Analysis 

Observed delays r depend on an orientation of the 
baseline vector and the vector of source position, as 
well as propagation effects and clock functions. The 
differences between observed delay r and predicted 
delay r c can be used for adjusting a vector of parame- 
ters Ap using the least squares: 



EOT 
— A PI =T -T C . (1) 

i 

The overview of the technique for computing pre- 
dicted delays can be found in Sovers, Fanselow, & Ja- 
cobs CS (1998). We used the software library VTD 4 
for computing the theoretical path delay. The ex- 
pression for time delay derived by Kopeikin & Schafer 
(1999) in the framework of general relativity was used. 
The displacements caused by the Earth's tides were 
computed using a rigorous algorithm of Petrov & Ma 

(2003) with a truncation at a level of 0.05 mm based 
on the numerical values of the generalized Love num- 
bers presented by Mathews (2001). The displacements 
caused by ocean loading were computed by convolving 
the Greens' functions with ocean tide models using the 
NLOADF algorithm of Agnew (1997). The GOT00 
model (Ray 1999) of diurnal and semi-diurnal ocean 
tides, the NA099 model (Matsumoto, Takanezawa, & 
Ooe 2000) of ocean zonal tides, the equilibrium model 
of the pole tide and the tide with period of 18.6 years 
were used. The atmospheric pressure loading was com- 
puted by convolving the Greens' functions with the 
output of the atmosphere NCEP Reanalysis numerical 
model (Kalnay et al. 1996). The algorithm of com- 
putations is described in full details in Petrov & Boy 

(2004) . 

The accuracy of the a priori model is not enough 
to resolve reliably multi-band group delays computed 
over the entire recorded bandwidth that have ambi- 
guities with spacing 4.5 ns. The ambiguity spacing 
of narrow-band group delays, that were computed for 
each frequency channel separately and then averaged, 
is 16 mks, but the precision or narrow-band delays is 
worse by the ratio of the width of the entire recorded 
band to the width of the individual channel, i.e. a 
factor of 14. The uncertainties of narrow-band group 
delays are in the range of 0.05-10.0 ns. Of those, 78% 
fall within less than 0.8 nsec, i.e. 1/6 of the multi-band 
group delay ambiguity spacings for 78% observations. 

4 Source code and detailed documentation is available at 
http : //astrogeo . org/vtd 
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Therefore, a direct substitution of narrow-band delays 
to multi-band group delays will leave wrong ambigui- 
ties for 1 /5 of the observations. The time- variable dif- 
ferences between narrow-band delays and multi-band 
delays due to instrumental effects will worsen the sit- 
uation even further. 

To circumvent the problem, we first made the so- 
lution using narrow-band group delays. The estimated 
parameters were coefficients of the 1st degree B-spline 
that model clock functions and atmosphere zenith path 
delay at each station, positions of the three new sta- 
tions, ATCA-104, MOPRA, CEDUNA, coordinates of five 
new sources, the pole coordinates and the UT1 angle. 
We formed new right-hand sides in a system of 

dr 



equations 1: 



-E 



r\ n 

OPi 



-Apr, where r 9 



are multi-band observed group delays and p™ are pa- 
rameter adjustments from the LSQ solution that used 
narrow-band delays. The accuracy of the a priori 
model with adjustments taken from the narrow-band 
delay solution turned out sufficient for reliable resolv- 
ing group delay ambiguities at all baselines. The ex- 
ample of the differences between observed multi-band 
group delays and predicted are presented in Figures 2- 
3. 




5 10 

Time in hours 

Figure 2 The differences between observed mufti-band 
group defays and theoreticaf delays for basefine ATCA- 
104/ CEDUNA predicted on the basis of the a priori model 
updated with results of the LSQ sofution that used narrow- 
band delays. 



During the preliminary fringe test, we used fre- 
quencies 23472, 23488, 23600, and 23616 MHz. This 
frequency setup resulted in ambiguity spacings 7.8125 ns. 
Analysis of the fringe test showed that multi-band de- 
lay ambiguities can be reliably resolved even with twice 
smaller ambiguities. The fringe amplitudes were very 
weak at baselines with HOBART26 and no fringes was 
detected at baselines with CEDUNA. We found that fre- 
quencies 23600-23632 MHz were near the edge of the 
receiver filter at HOBART26 and beyond the filter at 
CEDUNA. We have adjusted the frequency setup at the 
experiment on 2007 June 24. The fringe test allowed us 
to determine positions of ATCA-104 and MOPRA with 
l-<7 errors ~ 1.6 m for the vertical component and 
~0.2 m for the horizontal component. We found that 
the a priori coordinate of ATCA-104 had an error of 




£ 1 



0) 
Q 



— 1 



5 10 

Time in hours 

Figure 3 The differences between observed multi- 
band group delays and theoretical delays for baseline 
MOPRA/HOBART26 predicted on the basis of the a priori 
model updated with results of the LSQ solution that used 
narrow-band delays. 



80 m, apparently due to a confusion of different pads 
used by the ATCA array. 

After successful resolution of group delay ambigu- 
ities the outliers were eliminated, and the additive 
baseline-dependent corrections to weights was eval- 
uated in a trial solution to make the ratio of the 
weighted sum of squares of residual to their mathemat- 
ical expectation to be close to unity. The RCP and 
LCP data were treated as independent experiments. 
All analysis procedures, group delay ambiguities reso- 
lution, outliers elimination and reweighting were made 
independently for the LCP and RCP datasets. 

Finally, the RCP and LCP data from the 12 hour 
run were used in a global least square (LSQ) solution 
together with all available VLBI data, 4431 twenty four 
hour experiments from 1980 through 2008. 

The estimated parameters belong to one of the 
three groups: 

- global (over the entire data set): coordinates of 
3094 sources, including 5 new sources observed 
during the experiment on 2007 June 24, posi- 
tions of 152 stations, including positions of all 
stations observed on 2007 June 24. 

- local (over each session): tilts of the local sym- 
metry axis of the atmosphere for all stations and 
their rates, station-dependent clock functions 
modeled by second order polynomials, baseline- 
dependent clock offsets. 

- segmented (over 0.33-1.0 hours): coefficients of 
linear spline that models atmospheric path de- 
lay (0.33 hour segment) and clock function (1 
hour segment) for each station. The estimates of 
clock function absorb uncalibrated instrumental 
delays in the data acquisition system. 

The rate of change for the atmospheric path de- 
lay and clock function between two adjacent segments 
was constrained to zero with weights reciprocal to 
1.1 ■ 10~ 14 and 2 • 10~ 14 , respectively, in order to sta- 
bilize solutions. The weights of observables were com- 
puted as w — 1/ y/a% + r f , where a is the formal 
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uncertainty of group delay estimates and rt is the ad- 
ditive baseline-dependent reweighting parameter. 

Since the wave propagation equations are differen- 
tial equations, their solution does not determine object 
coordinates and their derivatives uniquely. As a result, 
the matrix of normal equations that emerges in solving 
equations 1 has incomplete rank. In order to overcome 
the rank deficiency, we impose constraints that require 
the net-rotation and net-translation of the estimates of 
positions and velocities of 35 selected sites with respect 
to positions from the ITRF2005 catalogue (Altimimi 
et al. 2007) to be zero, and the net-rotation of 212 
sources with respect to the ICRF catalogue (Ma et al. 
1998) be zero as well. 

Results of the solution are shown in Tables 3 and 
4 5 . The column (u) shows the number of observations 
used in the solution. The column (s) shows the num- 
ber of scheduled observations. The 1-a uncertainty 
of the vertical position of ATCA-104 and MOPRA are 
12 mm, and the uncertainty of their horizontal po- 
sition is 4 mm. The 1-a uncertainty of the vertical 
position of CEDUNA is 28 mm, and the uncertainty of 
its horizontal position is 6 mm. 

The reference position for the ATCA is set at ob- 
serve time to one of the 45 antenna pad positions. 
When the ATCA observes in the phased array mode, 
the signal from the reference antennas is delayed to a 
constant value depending on configuration, typically 
about 50 mks. This signal is mixed with the signals 
from other antennas of the array delayed in such a way 
that they alway arrive in phase with the signal from 
the reference antenna. This variable delay is computed 
in real time on a basis of the geometric model that in- 
corporates relative positions of antennas in the array, 
coordinates of the observed sources and time. The off- 
set between pad positions is known to millimeter preci- 
sion from local surveying and at any time the antennas 
are positioned on the pads to better than 10 mm pre- 
cision, so the position derived here for the ATCA can 
be converted to any ATCA reference pad position. Ta- 
ble 2 lists the nominal offsets for each ATCA pad with 
respect to the CATW104 pad. 

All five new sources were detected and their coor- 
dinates were determined with the 1-a semi-major axis 
of the error ellipse 2-5 mas according to Table 4. The 
last column "s" shows the number of scheduled obser- 
vations, the column "u" shows the number of obser- 
vations used in the solution, and the column: corr" 
shows correlation between estimates of right ascension 
and declination. 

4.2 Correlated Flux Density Esti- 
mates 

Amplitudes of the cross-correlation function are used 
for estimation of the correlated flux density of observed 
sources. System temperature was measured at all sta- 
tions during the experiment, either before each scan 
or continuously. The correlated flux density F corr is 



Table 2 Nominal offsets of ATCA pads with respect to 
the CATW104. These offsets were determined by the geo- 
detic survey when the ATCA was commissioned and they 
presented here for the reference purpose only. 

Pad name East offset North offset 
in meters in meters 



f~i A mi T 7"0 

CATWU 


1 COI QQ7 


U.UUU 


CATW2 


lool.zzo 


A AAA 
U.UUU 


CATW4 


1 CQA CI Q 

looU.Olo 


A AAA 
U.UUU 


CATWO 


i cnn aaa 
loUU.UUU 


A AAA 
U.UUU 


CATW8 


44oy.o88 


A AAA 
U.UUU 


y~-> A mi T T 1 (\ 

CATW1U 


1 A QQ 77 
14oo. I/O 


A AAA 
U.UUU 


CATW12 


1 A AQ 1 d A 

14U8.104 


A AAA 
U.UUU 


CATW14 


16 ( Cool 


A AAA 
U.UUU 


y~i A mi T T 1 C\ 

CATW lo 


io4o.yoy 


A AAA 
U.UUU 


C ATW3 2 


1 1 no a A 1 
11UZ.U41 


A AAA 
U.UUU 


CATW45 


nno A£i 

yuo.uoi 


A AAA 
U.UUU 


CATW64 


0lz.z40 


A AAA 
U.UUU 


CATW84 


onfi 1 OO 

oUo.lzo 


A AAA 
U.UUU 


CATW98 


yi.oo < 


A AAA 
U.UUU 


CATW 100 


Ol.ZZO 


A AAA 
U.UUU 


CATW102 


on fii o 

oU.Olo 


A AAA 
U.UUU 


LAI WTU4 


U.UUU 


n nnn 
U.UUU 


CATW 106 


— oU.Olz 


A AAA 
U.UUU 


CATW 109 


no con 


A AAA 
U.UUU 


CATW 110 


— yi.ooo 


A AAA 
U.UUU 


CATW 111 


1 fl7 1 A Q 
— 1U( .116 


A AAA 
U.UUU 


CATW112 


lOO A A n 

— lzz.44y 


A AAA 
U.UUU 


CATWl 13 


1 07 7CC 

— lo I . too 


A AAA 
U.UUU 


CATW124 


— oUo.lzz 


A AAA 
U.UUU 


CATW125 


Q01 /IOC 

— ozT.4Zc> 


A AAA 
U.UUU 


CATW128 


Qfi7 Q A 7 

—6vl .611 


A AAA 
U.UUU 


s~, a mi ir1 OH 

CATWlzy 


OQO RKQ 
— OOZ.DOO 


A AAA 
U.UUU 


CATW 140 


rri non 

— 551. OzO 


A AAA 
U.UUU 


CATW147 


— OOO.lDO 


A AAA 
U.UUU 


CATW148 


£J70 Add 


A AAA 
U.UUU 


CATW163 


ono n£i 
— »Uo. Ubl 


A AAA 
U.UUU 


CATW168 


070 cno 

— y < y.oyz 


A AAA 
U.UUU 


CATW172 


— 1040.816 


0.000 


CATW173 


-1056.122 


0.000 


CATW182 


-1193.877 


0.000 


CATW 189 


-1301.020 


0.000 


CATW 190 


-1316.326 


0.000 


CATWl 95 


-1392.857 


0.000 


CATW 196 


-1408.163 


0.000 


CATW392 


-4408.163 


0.000 


CATN2 


-30.612 


30.612 


CATN5 


-30.612 


76.531 


CATN7 


-30.612 


107.143 


CATNll 


-30.612 


168.367 


C ATN 1 4 


-30.612 


214.286 



5 A machine readable version of this table is also acces- 
sible in the electronic supplement. 
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Table 3 Position in meters of the three new stations on epoch 2007 July 24. 



Station 


X 


Y 


Z 


ATCA-104 


-4751640. 182 ±0.008 


2791700.322 ±0.006 


-3200490.668 ±0.007 


CEDUNA 


-3753443.168 ±0.017 


3912709.794 ±0.017 


-3348067.060 ±0.016 


MOPRA 


-4682769.444 ±0.009 


2802618.963 ±0.006 


-3291758. 864 ±0.008 



Table 4 Coordinates of five new sources at the J2000.0 epoch. 



Source 






a 




5 


(corr) 


(u) 


(s) 


0100-760 


0l h 


02 m 


18 s .6609 ± S .0018 


-75° 46' 


51 ".730 ± 0".003 


-0.47 


24 


32 


0219-637 


02 h 


20 m 


54 s .1722 ± 010004 


-63° 30' 


19".387±0".002 


-0.23 


31 


32 


0333-729 


03 h 


32 m 


43 S .0009±0 S .0011 


-72° 49' 


04". 521 ± 0".004 


-0.30 


18 


32 


1941-554 


19 h 


45 m 


24 s .2477 ± S .0007 


-55° 20' 


48".949 ± 0".006 


-0.30 


5 


16 


2140-781 


2l h 


46 m 


30 s .0694 ± S .0003 


-77° 55' 


54". 735 ± 0".002 


0.80 


25 


25 



proportional to the amplitude of the cross correlation 
A 



9i 92 A c 



(2) 



where T al and T s * 2 are system temperatures at sta- 
tions 1 and stations 2 corrected for the extinction in 
the atmosphere, Di and D 2 are the a priori elevation- 
dependent antenna gains for both antennas of the base- 
line, and gi , g 2 are empirical a posteriori multiplicative 
gain corrections. 

The a priori gains were obtained from prior ded- 
icated observations of bright sources. Its dependence 
on elevation was modeled with a polynomial. 

Recorded system temperature was considered as a 
sum of two terms: the receiver temperature T rec and 
the contribution of the atmosphere: 



■ + T a : 



i( e )l 



(3) 



where T a tm is the average temperature of the atmo- 
sphere, a is the atmosphere opacity, and m(e) is the 
wet mapping function: the ratio of the neutral at- 
mosphere non-hydrostatic path delay at an elevation 
e to the atmosphere non-hydrostatic path delay in 
the zenith direction. We omitted in expression 3 the 
ground spillover term that was not determined for 
these antennas. We set T atm to 280K, and evaluated 
the receiver temperature and the opacity factors for 
each antenna by fitting them into records of system 
temperatures with the use of the non-linear LSQ. The 
results are presented in Table 5. Column 2 contains 
the averaged adjusted System Equivalent Flux Den- 
sity (SEFD) at elevation angles > 60°. Column 3 
contains the multiplicative error factor of the SEFD 
estimate. Column 4 and 5 contain the estimates of 
the receiver temperature and the atmosphere opacity 
in the zenith direction. The system temperature di- 
vided by e~ am( - e \ T* , is free from absorption in the 
atmosphere. 



Table 5 Parameters of the radio telescopes determined 
from the experiment. 



Station 


SEFD Jy 


m.e.f. 


T rec K 


opacity 


ATCA-104 


820 


0.07 


25 


0.17 


CEDUNA 


9100 


0.09 


164 


0.10 


HARTRAO 


8200 


0.18 


206 


0.06 


HOBART26 


1800 


0.08 


144 


0.50 


MOPRA 


780 


0.07 


26 


0.12 


PARKES 


890 


0.10 


60 


0.06 


SESHAN25 


11000 


0.10 


119 


0.27 



The a priori antenna gain D and/or T ays may have 
a multiplicative error. The corrections to antenna gain 
were evaluated by fitting the correlated amplitude to 
the flux density of sources with known brightness dis- 
tribution. Of 97 sources detected in our experiment, 
38 objects were imaged in the VLBA K-band exper- 
iment on 2006 July 09 under the K/Q band survey 
program 6 . The brightness distributions of these ob- 
jects and many others were made publicly available 
by A. Fey 7 . We have computed the logarithms of the 
predicted correlated flux density at moments of obser- 
vations for these 38 sources considered as amplitude 
calibrators. We used them for determining logarithms 
of corrections to gains by the iterative LSQ procedure 
to fit to the logarithms of observed correlated ampli- 
tudes according to equation 2. The iterative procedure 
rejected 9 calibrators since the ratio of the measured 
correlated flux density to the correlated flux density 
from the map either exceeded 1.66 or was less than 0.6. 
The rms of the diffrences of these ratios from 1.0 over 
remaining 29 objects was 0.30. Although the corre- 
lated flux density of the calibrators may change for one 
year between epochs of observations due to image evo- 
lution, the large number of calibrators provides rather 
robust estimates of gain corrections. The estimates 
of the mean system equivalent flux density, defined as 



6 G.E. Lanyi et al, submitted for publication to Astro- 
nomical Journal, 2009. 

7 Availablc at http://rorf.usno.navy.mil/rrfid.shtml 
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g T sya e" m ' e ' / D — the parameter that characterizes 
the sensitivity of a radio telescope, are presented in Ta- 
ble 5 for elevations [60°, 90°]. Since the corrections to 
gains are multiplicative, their errors are characterized 
by a multiplicative error factor (m.e.f.). 

Using gain corrections, we have computed the cali- 
brated flux density for other 68 detected sources. Since 
we have too few observations of each individual source, 
we made no attempt to produce images. We have 
computed the mean weighted correlated flux density 
in three ranges of the baseline projection lengths: 0-50 
megawavelengths, 50-150 megawavelengths, and 400- 
800 megawavelengths. The results are presented in 
Table 6. The errors of correlated flux density esti- 
mates for point-like sources are determined by errors 
of amplitude calibration that are the mean m.e.f. from 
Table 5, i.e. ~ 15%. 

5 Conclusions and Future Ob- 
servations 

The K-band geodetic VLBI experiment with using the 
LBA network and SESHAN25 turned out highly success- 
ful. Position of ATCA-104 and MOPRA were determined 
with l-o" formal uncertainty 12 mm for the vertical 
component and 4 mm for the horizontal components. 
The catalogue of site positions and velocities in our 
solution does not have net rotation and net transla- 
tion with respect to the IRTF2005 catalogue, there- 
fore our estimates of site positions are consistent with 
the ITRF2005. Since relative positions of 44 ATCA 
stations were previously measured with ground sur- 
vey, using our coordinates of ATCA-104 we have de- 
rived absolute positions of other 43 pads. The table of 
absolute positions of the ATCA antenna stations can 
be found at http://astrogeo.org/lcs/coord. Posi- 
tion of CEDUNA was determined with the uncertainty 
28 mm for the vertical and 6 mm for the horizontal 
components. Worse position accuracy of CEDUNA is 
explained by 10 times worse antenna sensitivity and 
the lack of the LCP data. 

We have demonstrated that VLBI experiments for 
absolute astrometry and geodesy at the array with a 
heterogeneous setup that uses the LBADR and Mark- 
5 data acquisition systems are feasible and can provide 
high quality results. 

We have demonstrated that the group delay am- 
biguities with spacings as small as 4.5 ns can be suc- 
cessfully resolved using the adjustments to the a priori 
model from the narrow-band delay solution. 

We have determined positions of 5 new sources 
never observed before with the VLBI with the 1-cr un- 
certainties of 2-5 mas. This result proves that the LBA 
can be used for absolute astrometry observations., 

Inspired by these astrometric results, we launched 
the X-band LBA Calibrator Survey observing cam- 
paign for determining positions and images of thou- 
sands of objects with declinations below —50°. The 
first observing session ran successfully in February 
2008. 

8 http : //astrogeo . org/lcs 



Table 6 The mean correlated flux density of observed 
sources in jansky, except 29 amplitude calibrators used for 
the gain adjustment, at three ranges of the baseline projec- 
tion length. 

Source # Obs Correlated flux density in Jy 
0-50 50-150 400-800 
MA MA MA 



0007+106 


10 


0.43 


0.48 




0047—579 


23 


1.65 


1.22 




0048—097 


13 


1.21 


1.31 




0048—427 


8 


0.57 


0.36 




0100—760 


g 


0.36 


0.18 




0104—408 


10 


3.46 


2.07 




0107—610 


5 


0.37 


0.37 




0131—522 


17 


1.01 


0.65 




0219—637 


19 


0.36 


0.31 




0230 — 790 


10 


0.68 


0.46 




0234+285 


16 


3.86 


3.15 


1.10 


0235+1 64 


15 


2.85 


3.60 


2.56 


0252 — 549 


g 


1.41 


0.73 




0302 — 623 


17 


1.74 


0.91 




0306+102 


23 


0.94 


0.99 


0.71 


0322+222 


g 


0.41 


0.63 




0332—403 


3 


0.53 


0.39 




0333 — 729 


11 


0.25 


0.24 




0333+321 


3 


1.03 




0.43 


0336 — 019 


27 


2.00 


1.49 


0.70 


0358+040 


17 


0.51 


0.51 


0.45 


0402 — 362 


13 


1.12 


0.67 




0438—436 


11 


0.67 


0.48 




0446+1 1 2 


4 


0.55 


0.96 




051 6 — 621 


g 


1.15 


0.97 




0537—441 


6 


8.58 


6.37 




0552+398 


3 


2.39 




0.89 


0727—115 


1 


2.20 






0736+017 


3 


1.60 


3.63 




0738 — 674 


1 


0.22 






1057—797 


29 


2.58 


2.03 


0.86 


1144+402 


2 






1.24 


1144—379 


3 




2.37 


1.13 


1156+295 


2 






0.84 


1324+224 


2 






1.31 


1334—127 


16 


7.13 


10.78 


5.57 


1349—439 


3 


0.40 


0.27 




1406—076 


10 


1.00 


0.95 


0.81 


1424—418 


10 


2.33 


2.04 


0.86 


1 51 1 —100 

-LxJJ-J. IvU 


27 


1.32 


1.10 


1.38 


1610—771 


15 


1.65 


0.72 




1619—680 


14 


0.74 


0.42 




1718—649 


7 


2.36 


1.09 




1732+389 


2 


0.66 




1.08 


1741 —038 


14 


3.99 


3.15 


2.49 


1 74Q+096 


16 


6.03 


8.41 


4.41 


1758 — 651 


18 


0.90 


0.72 


0.43 


1 800+440 


3 


1.47 




0.61 


1815 — 553 


11 


0.73 


0.45 




1831—711 


9 


1.62 


0.99 


0.46 


1925—610 


16 


0.57 


0.30 




1941—554 


4 


' ' * 


0.24 




i nc/i 100 
1954— 3oo 


16 


2.37 


1.37 




1958-179 


25 


2.44 


1.79 


1.53 


2030-689 


24 


0.52 


0.30 




2052-474 


10 


2.58 


2.03 


0.90 


2059-786 


10 


1.03 


0.65 




2121+053 


18 


1.16 


0.93 


0.50 


2131-021 


23 


2.02 


1.73 


1.53 


2140-781 


15 


1.15 


1.00 


0.97 


2142-758 


12 


0.40 


0.26 




2145+067 


26 


5.93 


5.21 


2.09 


2227-088 


16 


4.45 


3.63 


4.37 


2236-572 


38 


0.80 


0.71 




2326-477 


36 


1.13 


0.80 




2344-514 


22 


0.27 


0.23 




2353-686 


10 


0.93 


0.65 




2355-534 


10 


1.96 


1.57 
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